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ABSTRACT 

Low temperature (-150°C) EPR spectra of the nitroxide generated 
photolytically in a single crystal of N-nitromorpholine are complicated 
by the presence of two magnetically distinct orientations of the N-0 
bond within the crystal. As the temperature is raised, the spectra are 
further complicated by motional effects. Since at -150° the molecule 
is essentially frozen in the crystal lattice with respect to the time of 
the EPR interaction, and since the nitrogen coupling tensor is markedly 
anisotropic, the observed EPR spectrum varies widely as the orientation 
of the applied field changes. The anisotropy of the nitrogen coupling 
tensor made it possible to unequivocally determine the orientation of 
the crystal with respect to the Poole magnetic field. A series of 
experimental spectra was then obtained for several combinations of 
crystal orientation and temperature. In addition, computer programs 
were developed which simulate low, intermediate, and high temperature 
EPR spectra. By using a combination of experimental spectra and these .- 
programs, it was possible to determine the angle between the two 


orientations of the N-0O bond in the crystal and also to determine the 


energy barrier to exchange. 
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ia INTRO BUe TION 


A. BACKGROUND 

EPR studies of paramagnetic species trapped in solids may be 
conveniently classified according to the matrix in which the species is 
caged. Some common classifications are powders, glasses (frozen 
solutions), host crystals and single crystals. Each class has its own 
advantages and disadvantages, depending on the particular experimental 
situation. In general, however, if the experimental situation lends 
itself to the use of single crystals, then it will be found that single 
crystal EPR studies are to be preferred over the other types. In particular, 
since the species to be studied is regularly oriented within the crystal 
matrix, it is possible to determine the directional dependence and 
magnitude of anisotropic effects such as the dipole-dipole interaction. 
A set of orthogonal crystal axes can be assigned on the basis of the 
outward appearance of the crystal, and then these crystal axes may be 
used to determine the orientation of the crystal within an applied magnetic 
field. In addition, direction cosine matrices can be experimentally 
determined which relate the principal axes of the coupling tensors to 


the crystal axes and thus to the direction of the applied field. 


EB PREVIOUS REonawk © oH 
Single crystal EPR studies of the nitroxide of N-nitromorpholine 


(hereafter referred to as morpholine N-oxide) by Bodnar [1] showed the 


10 








presence of two magnetically distinct sites at room temperature. The 
sites were characterized as a Single nitrogen nucleus and two equivalent 
pairs of protons. Interchange between the axial and equatorial protons 
at room temperature was sufficiently rapid to cause an averaging effect 
so that hyperfine coupling for each pair of B protons could be labeled 
as A, = 5 a 

Bodnar developed the following set of hyperfine coupling tensors 
for morpholine N-oxide (expressed as column vectors): 

Table I: Hyperfine Coupling Tensors 
‘ for Morpholine N-oxide 


Nucleus Principal Values and Direction eosines 





AWe= 8 .ore A. = Seer eee! 3 GS 
yY 22 


XX 
x 0.83 HON 0.54 
i y 0.23 0.97 ae 
ee =(0)).52 0.15 0.84 

Awe 13.6 fh = ee Cc te ene 
eK VY ZZ 
x 0.84 0.42 0.34 
4H y 0.01 0.62 -0.78 
z  -0.54 0.66 0.52 


Windle, Kuhnle, and Beck [2] studied the EPR spectra of morpho- 
line N-oxide in methylene chloride over a temperature range of -130° to 
+ 120°C and found that the sum of the axial and equatorial hydrogen 


coupling constants increased with decreasing temperature. 


ie 








ic. OBJ ECH @F THissREoEeR@ rH 

The primary goal of this research was the development of tools 
which would simulate the behavior exhibited by complex paramagnetic 
systems of the sort typified by morpholine N-oxide. These tools could 
then be applied to the secondary research aims of: 

ils Refinement of the observed coupling constants of the 
morpholine N-oxide radical, including variations with temperature 
and crystal orientation. 

oe Simulation of EPR spectra in the absence (or presence) of 
motional effects. 

oo Determination of the angle between the magnetically distinct 
Orientations of the N-0O bond in morpholine N-oxide at low temperatures. 

4. Determination of the energy barrier to nuclear exchange for 
morpholine N-oxide in a matrix of N-nitromorpholine 

a: Identification of the directions and magnitudes of the dipole- 
dipole and Fermi Contact interactions which occur in complex para- 


magnetic systems of this sort. 
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Ils THEORY AND DEVELOPMENT eRe trea holes 


A. THE APPROPRIATE HAMILTONIAN 

Carrington and McLachlan [3] characterize the spin Hamiltonian 
for the EPR interaction of one unpaired electron near a single nucleus 
as: 

+ Ss ee an 
oie we ecm ee chee a 
The terms in the Hamiltonian respectively describe the contributions 
to the total energy which are due to interactions between the applied 
field and the unpaired electron, the electron and the nucleus, and the 
applied field and the nucleus. The first term represents the interaction 
st a 
between the applied field H and the spin moment S of the unpaired 
electron (electronic Zeeman energy). [2 is the Bohr magneton, and G is 
the spectroscopic splitting tensor. The second term is a combination 
of two terms and represents the hyperfine interaction between S and the 
nuclear spin ig (Fermi Contact and dipole-dipole interactions.) i is 
the hyperfine interaction tensor. The third term represents the inter- 
action between the applied field and the nuclear spin (nuclear Zeeman 
energy). em is the spectroscopic splitting factor for the nucleus, and 
3 F is the nuclear magneton. 
i The g-Tensor and the Electronic Zeeman Energy 
In the absence of spin-orbit coupling, the g-tensor is 


we 
isotropic and S represents the true spin of the electron. However, 


es 








when the electron possesses both spin and orbital angular momentum, 
the g-tensor is anisotropic, and S no longer represents the true electron 
spin since the spin is not oriented precisely along the field direction. 
An acceptable way to treat this problem is [4] to assume the existence 
of an "effective" field, such that the effective field direction is along 
the direction of quantization of the electron spin wector. The development 
is as follows: 
The principal values of the g-tensor can be determined by 


choosing as system axes the principal axes of the tensor. Ifg , oe 
soe 


g_, are the principal values of the g-tensor, and oe y o, are the 


direction cosines of the applied magnetic field in the principal axis 


“n 


oad = ; mN 
system, thenH:> g=H(g y t+ g See Ie Sa kk). Thesvecion 


Galea meas 


> «=n 
quantity in parentheses can be more simply expressed as op h, where 


ff 
epee ee as oa | 
[4] oa lig, ¥ .) fs g,. vy) ne (g. ¥.) weaned ao bie unat 


Ee EX 
vector along which the electron spin is quantized. ThenH * g= 


fae h=g Hi . Finally, the electronic Zeeman interaction 
eff Stic 
becomes 2 = © oe, Hoes amen 
2 The A-Tensor, the Fermi Contact Interaction and the 


Dipole-Dipole Interaction 


The second term of the spin Hamiltonian may be separated 

mato an isotropic part aS « I and an anisotropic part S- T- I. Thus 
a 3 
the A-tensor can be represented as: A=T tal. 


The isotropic part aS - I arises from the Fermi Contact 


interaction which represents the hyperfine coupling between the nuclear 
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spin moment and the magnetic field produced at the nucleus by the 
spinning electron. The contact interaction can only occur when the 
electron has a finite probability censityeat Encummeloucmiaye 
as cad > 
The anisotropic part S - T - Iis due to magnetic coupling 
between the magnetic moments of the electron and nucleus. This 


dipolar portion of the Hamiltonian is given by 


—e 
ee 


—_ oie ; aE eee _ oS ~ = 
ee ee 9 BLoP es 3(S uo | 
r if 


ob 
where r is the radius vector between the two moments, andr is the 
distance between them. The dipolar portion of the Hamiltonian averages 


out to zero whenever the electron cloud is spherical [3]. 


oe Combined Effects of the Fermi Contact, Dipole-Dipole and 


Nuclear Zeeman Interactions 
The second and third terms of the spin Hamiltonian may be 
combined as 


—> 


—- — ad = . 
nie 8 een ee 


> Lb 


~~ 
S e 
> wy 
In addition, S may be represented as mh, where mo is the spin quantum 
—s 
number of the electron (+3) and h is a unit vector in the direction of 


the applied field. Substituting above and rearranging gives 





ee me eg ee 
n By n ae g B n Pn eff 
n n 
4, Resultant Hamiltonian for the General Case 


Using the expressions developed in sections 1] and 3 above, 


and expanding to the case of a single electron interacting with several 


LS 








nuclei, the spin Hamiltonian becomes 
a > —> _ > es P —- 
1 eee Svate : ; iS. 2 aoe I} 
where I. is the spin of the i-th nucleus. 
i 


. PROGRAM ESRPLT 

The Hamiltonian given in section A-4 above was used as the basis 
for a computer program which would predict and plot complex EPR spectra 
in the absence of motional effects. The mechanics of this program, 


called ESRPLI, are given in Appendix I. 


c.. PROGRAM NITROX 

A second computer program was developed which was based ona 
simpler treatment of the EPR interaction than the approach used in ESRPLT. 
This program, called NITROX, makes use of the high field approximation 
[5, 6], which assumes that the nuclear Zeeman effects (NMR transitions) 
are insignificant compared to the electronic Zeeman effects (EPR 
transitions). 

NITROX was designed to compute coupling constants for a nitroxide 
as a function of the direction of the applied field. Coupling constants 
were formulated using a variation on the method employed by Lefebvre 


and Maruani [6] 


where A is the computed coupling constant, A ,A_,A__ are the 
Moc ty eZ 
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principal values of the input coupling tensor, and s. O and io 
are the direction cosines of the applied field in the principal axis 
system. (Note the similarity of this approach to the method used to 


.) 


compute Daeg 
NITROX also computes the energies and intensities of the lines 
of the first order EPR spectrum, plots the spectrum, and punches cards 
which contain energy-intensity data in a format which is required by 
program GEX. This program is adequate for the simulation of EPR 


spectra in the presence or absence of motional effects. The mechanics 


of program NITROX are given in Appendix II. 


ie) PROGRAM GEX 
_A third computer program, developed by Tolles [7], simulates the 


shape of an EPR spectrum as a function of the rate of exchange. 
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III. EXPERIMENTAL 


A. IDENTIFICATION OF CRYSTAL AND PRINCIPAL AXES 
Ne Assignment of Crystal Axes 

The crystal axes assigned for crystals of N-nitromorpholine 
were those used by Bodnar [1]. The orientation of these crystal axes 
relative to the physical appearance of the crystal is depicted in 
Figure l. 

ae Crystal Orientation in the EPR Spectrometer 

The method used to suspend the crystal between the pole 
pieces of the EPR spectrometer is shown in Figures 2 and 3. A very 
thin layer of Kel F grease was used to make the crystal adhere to the 
flat ground surface of a quartz rod. The intent was to have the crystal 
aligned so that the + Z axis pointed toward the observer, and the + Y 
axis pointed downward when the quartz rod was viewed in a direction 
perpendicular to the flat ground face with the long extension of the rod 
pointed upward. Adjustments in the angle @ were accomplished by 
using a microscope equipped with a 360 - degree turntable and a cross- 
hair eyepiece. As the crystal was viewed along the + Z axis, clockwise 
rotation of the crystal about the Z axis corresponded to a counter- 
clockwise movement of the applied field, increasing @. When the 
crystal was suspended between the pole pieces of the spectrometer, a 


pointer attached to an extension of the quartz rod indicated 8. Variation 
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Figure 1. Crystal Form and Assigned Axes 
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Figure 2. Crystal Orientation on the Quartz Rod 





Figure 3. Alignment of the Crystal in the 
Maegenetic Field of the EPR Spectrometer 


20 








of 8 could be accomplished by moving this pointer. When viewed from 
the top, (along the -y axis), counter-clockwise rotation corresponded 
to a + change in 8. 

It can readily be seen from Figure 1 that some ambiguity could 
easily arise in the determination of crystal orientation relative to an 
applied magnetic field. Once crystal axes have been assigned, it is 
theoretically possible to ascertain crystal orientation unequivocally. 
However, it is often difficult to discern from the physical appearance 
of the crystal whether or not one of the axes has been inverted. This 
ambiguity can be readily resolved, however, by making use of previously 
tabulated values of coupling constants [1] and by knowing the possible 
confusing orientations of the crystal. Some selected values for the 
coupling constants of morpholine N-oxide at room temperature are 
tabulated in Table II versus the orientation of the applied field (polar 


coordinates) relative to the crystal axes. 


Table II: Room Temperature Coupling Constants 
for Morpholine N-oxide 


Be) i) ane (Gauss) a (Gauss) 
0 -45 Z7.4 120 
30 0 cee’ 11.4 
60 | -45 a) ey 
120 0) Bi Loeee 
120 -45 ees Pa 


Za 








For any given orientation of the crystal, the possible axis 


transformations are 


1. A 180-degree rotation about the X axis; 
Zs A 180-degree rotation about the Y axis; 
Sey A 180-degree rotation about the Z axis. 


Since morpholine N-oxide is symmetric with respect to a reversal 
in the direction of the applied field, there is a pair of equivalent polar 
coordinate angles corresponding to a 180° rotation about any one of 
the crystal axes. These pairs of polar coordinate angles which arise 
from axis inversion can be seen to stem from changes in the signs of the 
elements of an identity matrix as the matrix is subjected to the operations 


stated above. Consider the identity matrix 


Oe 
O.4i 6 (0 
One awe 


where the first, second, and third row non-zero elements correspond 


respectively to 


x = sin 8 cos @ 
y = sin 8 sin @ 
Z= cos 8 


It can now be seen that the following relationships hold 


i For a 180° rotation about the x axis we have 
_H - -H 
] 8@—>180° - 86 -] . G—>Q 
-] Or l 
~]) G—- - @ 1 G6—»180° - 7 


ZZ 








lie For a 180° rotation about the y axis we have 


z H . -H 
-] G@—»-180° - 86 1 B—» 6 
] or -] 
-l} g—s180° - Gg 1 g— - g 
Ce For a 180° rotation about the z axis we have 
eee ia = ll 
-] G—> 9 1 G@—>180° - 86 
-] or ] 
] g—>180°+9 -] g— g 


In general, the particular set of 8,% selected to represent a rotation 
about a given crystal axis was chosen such that it was physically 
possible to accomplish the orientation change using the standing set-up 
of the EPR spectrometer. 
Ore Principal Axis Systems 
The directions, in the crystal axis systems, of the principal 
axes of the several coupling tensors can be extracted from the appropriate 
matrix of direction cosines. The directions of the principal axes of the 
nitrogen coupling tensor were determined uSing the matrix of direction 
cosines found on page l]l. The orientations of these principal axes 
relative to the crystal axes are shown in Figure 4. 
4. The Angle Alpha 
The angle alpha was defined as one-half of the angular 
displacement which exists between the orientations of the N-O bond. 
This assignment was a matter of convenience since the basic direction 
cosine matrix for the nitrogen coupling tensor was the result of room 


temperature observations. At room temperature the exchange rate is 
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high enough so that what is seen is the averaged orientation of the 
principal axes (i.e., alpha = 0 on the average). When the exchange 
rate is slow relative to the time of the EPR interaction, each orientation 
of the N-0O bond results in a separate contribution to the observed 
spectrum. Thus alpha was taken as the angle between the "average" 
Orientation at room temperature and either of the two contributing 


Orientations at low temperatures. 


Ie. EXPERIMENTAL SPECTRA 


ile The EPR Spectrometer 


The EPR spectrometer used was a Varian model V-4502-13 
(X-band), with Fieldial V-F-2503 magnetic field regulation for the nine- 
inch magnet. The modulation frequency used was 100 KC. The EPR 
cavity was contained within a Dewar flask to facilitate low temperature 
studies. Temperature regulation was accomplished by directing a stream 
of cold nitrogen gas through a thermostatted heating unit at the base of 
the Dewar flask and into the cavity. Cavity temperature was measured 
by means of a copper-Constantan thermocouple. 

2 Synthesis of Morpholine N-oxide 

The morpholine N-oxide radical was generated by irradiating 

Single crystals of N-nitromorpholine with ultraviolet light. The source 


used was a Christie Xenolite Model UF1OJKK xenon-mercury vapor lamp. 


Zo 


oe Spectra Recorded 
Spectra were recorded for six crystal orientations at 

temperatures ranging from + 30° to -150°C. The orientations chosen 
were defined by the room temperature principal axis system of the 
nitrogen coupling tensor. The magnetic field was oriented successively 
along the Z, Y and X principal axes and the XZ, XY, and YZ diagonals. 
The relationships of these principal axes to molecular structure was 
proposed by Bodnar [1]. The Z axis is directed along the nitrogen 2P+ 
orbital where the unpaired electron is expected to be most highly 
localized. The Y axis is directed along the N-O bond, and the X axis 
is perpendicular to the plane containing the N-O bond and the nitrogen 
2P<, orbital. These orientations were chosen because it was expected 
that they would be most useful in the determination of the angle alpha. 
The temperature range was dictated by a combination of crystal melting 
point (51° - 52°C), the temperature regulation system, and the observable 
changes in the EPR spectra as temperature was lowered. Table III gives 


a summary of spectra recorded. 
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Table III: Summary of Experimental Spectra 


Magnetic Field Orientation Temperature 
Crystal Axes (8,2) niga Principal Axes ae 
32.00, ono) Z axis Wa, 0, -~1o5j=—Zo., 00 alow 
(S237 99.90) Y axis 20,0,-50,-100,-150 
15; 15.3) X axis Oeaue -o, 20 > Operate Oe 
W275, 8,5) XZ diagonal 29,20,0,-15,-50,150 
17.5, =-37 20) XY diagonal eo 10, -15,=50 ,-150 
(60.5, -54.5) YZ diagonal Sere oO oe) 
-50,-150 


ic SIMULATED SPECTRA 


ee The g-Tensor 


The principal axes of the g-tensor were assumed to coincide 
with the principal axes of the nitrogen coupling tensor. The principal 
values of g used were determined by Griffith, Cornell, and McConnell [8]. 
The values of et ae and ae used were 2.0061, 2.0089, and 2.0027 
respectively. (Note the interchange in assignments of the X and Y axes.) 

Le. Determination of the Angle Alpha 

The spectra which were simulated for the determination of 
alpha were those which were estimated to be the most likely to show a 
large variation in the simulated spectrum for a small change in alpha. In 
addition, it was assumed that any changes in molec len orientation 


could be represented as angular rotations about a vector perpendicular 
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to the plane containing the ie - 2P~ orbital and the N-O bond. 
Consequently, the magnetic field orientations chosen were (principal 
axis system): 

a. Along the Z axis (nitrogen 2P4y orbital), 

Di, Along the Y axis (N-0O bond), 

cn Along the YZ diagonal (midway between the N-0O bond 

and the 2P 4, orbital in the YZ plane). 
oe Determination of the Rate of Exchange 

An inspection of experimental spectra revealed that the 
series of spectra which showed the most clearly distinguishable changes 
with changing temperature were those for the case where the applied 
field was directed midway between the X principal axis and the 2P+y 
Sroital im the XZ plane (6 = 77.5°, @ =8.5° crystal axis system). As a 
result this orientation was chosen for exchange studies. Program NITROX 
was used to compute the line energies and intensities which make up the 
spectrum for the case of no exchange, and the punched output was used, 
in turn, as input data for program GEX. GEX was then used to simulate 


the spectrum changes as a function of the rate of exchange. 
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IV. RESULTS AND OBSERVATIONS 


a COUPLING CONSTANT VARIATION = EXPERIMENTAL SPECTRA 
The hyperfine coupling constants which were observed for 


morpholine N-oxide at -150°C are summarized in Table IV. 


Table IV: Hyperfine Coupling Constants for 
Morpholine N-oxide at -150° 


Nucleus Values Observed (Gauss) 
Ayy Avy ‘D2 


ea Tsay, (ees) 30n0 
cee 230 Ue 8 Zo 
equatorial H+ ; ions Dads Zee 


Note 1: Hydrogen coupling constants are those observed along 
An ee oe 
the N principal axes and are not true principal values. 
14 
ie Observed N Coupling Constants 
14 . 
The N coupling constants observed for magnetic field 
Orientations along the principal axes of the nitrogen coupling tensor 
are recorded as a function of temperature in Table V and are depicted 
ita@ee igure 5. 
Two interesting observations may be made regarding this 
data. First, the average of the principal values of the coupling constants 
at -150°C is 16.26 gauss, which is in good agreement with the isotropic 


value of 16.8 gauss (-100°C) reported by Windle, Kuhnle and Beck [2]. 
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Second, and perhaps even more interesting, is the observation that the 
average of the coupling constants at 20°C and -150°C are, for practical 


purposes, the same (16,27 vs 1le.26 Gauss). 


1 
Table V: N é Coupling Constants as a 
Function of Temperature 


C.C (gauss) TALC) 

Bees 20 

8.3 0 

Avy 7 ot -25 

PS -100 

Tele? -150 

OF) 20 

9.0 0 

A : Sd 

eo 

- -100 

7.8° -150 

3: lS 20 

52 ale 0 

Ago 325 -25 

3203 =50 

Soc -150 

Notes: 

i Some coupling constants were impossible to determine from 


experimental spectra because of poor resolution. 
ee . This coupling constant was determined using a combination of 


experimental and simulated spectra. 
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Dee Observed Hydrogen Coupling Constants 


14 
Particular observations, such as those made for the N 


coupling constants, cannot presently be made regarding the coupling 
constants for the B protons. A few general observations, however, 
are possible from the data which was recorded. 

a. Anisotropy in the coupling constant was definitely 
observable as the temperature was lowered. Although the separate 
principal axis systems for the axial and equatorial protons were not 
determined at low temperatures, spectra recorded for magnetic field 

. Ns cuca — 
orientations along the N principal axes indicate the presence of 
distinct coupling tensors for the two kinds of {3 protons. 

De The values of the axial proton coupling constants 

ee | | 
observed along the N principal axes did not vary widely (21.8 to 
2323 gauss) relative to the equatorial coupling constants for the same 
ee@entations (7.3 to 2:2 gauss). The lower figure for the equatorial 


proton coupling constant was estimated based on experimental and 


Simulated spectra. 


iB. SEPLC IRL M VARIATION Will TEMPERATURE 
As the temperature was lowered from 20° to -150°C, the recorded 
Spectra showed several distinct changes which were independent of the 


direction of the applied field. 


oi 








i Resolution became poorer as the temperature was lowered. 
Lines broadened and disappeared, or coalesced into new or already 
existing lines. 

Ze Coupling constants associated with N a increased or 
decreased, depending on the orientation of the applied field. This 
effect was noted previously and will not be discussed further here. 

The decrease in spectrum resolution which occurred as temperature 
was lowered below 20°C is directly related to the rate of interconversion 
between alternative conformations of the radical species. The effect 
on the spectrum of the rate of exchange of the axial and equatorial 
protons is shown by the correlation diagram in Figure 6. The base 
spectrum for a low rate of exchange is that which would result from the 
coupling of a single ee nucleus with two inequivalent pairs of protons 
(as is the case for morpholine N-oxide). However, because the 
existence of two separate orientations of the N-O bond is distinguishable 
to the EPR spectrometer when exchange is slow, the correlation diagram 
represents only part of the reason for poor resolution. The simulated 
spectra in Figure 7 show another complicating factor. The spectrum at 
the top of the page is the one which would exist if only a single 
Orientation of the N-O bond were distinguishable. The superimposed 
spectra at the bottom of the page show what actually exists. At low 
exchange rates the observed spectrum is the composite of the super- 
imposed spectra in Figure 7, and only nine distinguishable peaks 


represent a total of fifty-four contributing lines (see Figure 9). When 
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Figure 6. Correlation Diérram Showings the Effect on the 
Observed Spectrum of the Rate of Exchange of tne 
Axial and Equatorial Protons. Bottom: slow 
Sucuence Oo. ranrc OxXCNense; maiale: tae. 
mediate case. 











Top: Simulated Low Exchange Rate Svectrum for 
a Single Orientation of the N-O Bond 

Bottom: Simulated Low Exchange hate Syect run 
for the Case Where Two Orientations of the 


N=OuBonG fost. 
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a high rate of exchange exists, molecular motion results in an averaged 
picture so that the top section of the correlation diagram is a good 
representation of the observable spectrum. An idea of the complexity 
of the observed spectrum at intermediate rates of exchange can be had 
by mentally superimposing slightly displaced correlation diagrams on 
the lower spectrum in Figure 7. (The 1:1:1 peaks would coalesce for 


the high rate of exchange case.) 


ic. DETERMINATION OF ALPHA 

A combination of spectra recorded at -150°C and simulated 
spectra were used for the determination of alpha. As was noted 
previously, magnetic field orientations were chosen to lie in the plane 
containing the n'4 ZR OY Ditaleome ener Ni —Usbend- 

After experimentation with both ESRPLT and NITROX, NITROX was 
chosen for further work since both programs produced identical spectra 
for the case of no exchange, and NITROX took about one-twentieth as 


much computer time. In addition, NITROX is the more flexible program. 


Le re Coupling Constants - Simulated Spectra 


The principal axis system used to represent the The coupling 
tensor was the one determined at room temperature. The existence of 
two orientations of the N-0 bond was simulated by displacements (by 
the angle alpha) away from these axes. Since the displacements were 
assumed to occur in the YZ principal axis plane, only the Y and Z input 


principal values could affect the predicted nitrogen coupling constant in 


36 





the simulated spectra. Taking the case where the applied field is along 
the Z principal axis, for instance, the computer determined the lines 
and intensities of the spectrum which would result from the application 
of a magnetic field (+) alpha degrees away from the Z axis in the YZ 
plane. It then computed the lines and intensities which would result 
from the application of the field (-) alpha degrees away from the Z axis 
in the YZ plane (field displaced toward (-)Y axis). Finally, it combined 
the two sets of lines and intensities into one composite spectrum. 


| 
AG enon Ii were measured 


Now the observed values of Agog VY 


to be 33.8 and 7.8 gauss, respectively, at -150°C. These values 
include the effects caused by alpha. The trigonometric relationships 

14 
which describe the variation in observed N Zand Y coupling constants 


with alpha and with the input principal values Ao and Avy are: 


Nl 


2 
LL, 


A__ (observed) = (A 33.8 gauss 


dy 
(observed) = (A : eee +A ee ) 
Avy No ZZ; 


Bites ae ane ) 


NO]|— 


(oa 0als st. 


II 


The variation of Ano (input) and Ay (input) with alpha which occurs 
for correct predicted output values is plotted in Figure 8. 
Although there are any number of combinations of alpha, Dog 
(input) and Ay (input) which will result in a correct simulation of 
the observed spectrum for orientations of the magnetic field along the 
Y or Z principal axes, there is only one combination which will give 


the correct spectrum when the applied field is along the YZ diagonal 


ieee section lll, 6.4. fOr Claritieation Of this Orientation). 
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Figure 8. Variation of N 4 Coupling Constants with Alphe. 
These Coupling Constants are Correct Input 
Values for Computer Simulation. 
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ae Alpha 


By using input coupling constants from Figure 8 and assuming 
isotropic axial and equatorial proton coupling constants (as measured 
or estimated from experimental spectra), it was possible to determine 


a value of five degrees for alpha. 


i, SIMULATED AND EXPERIMENTAL SPECTRA 

Sets of experimental and simulated spectra for four different 
applied field orientations are compared in Figures 9, 10, ll, and 12. 
Regarding proton coupling constants, it can be seen that only in the 
case where the applied field was perpendicular to the 2P, orbital/N-0 
bond plane (Figure 12) was there a measurable (and large) equatorial 


proton coupling constant. 


es THE BARRIER EO EXGHANGE 

The change in observed spectra which can result from a relatively 
small change in temperature can readily be seen from Figure 13. The 
magnetic field orientation used to obtain the spectra in Figure 13 was 
the one used in programs NITROX and GEX to determine the energy barrier 
to exchange. Because of the general similarity of all of the spectra for 
any given orientation for temperatures below -50°C, the most reliable 
data was ana for temperatures greater than -50°C. The exchange 
rate data from program GEX were used to construct the Arrhenius rate 
plot depicted in Figure 14, where ~T is the rate of interconversion 


between conformers. The data and the slopes of the two lines shown 


So 





in Figure 14 were used to calculate the free energy, the enthalpy, and 

the entropy of activation for the process of interconversion. The appropriate 
formulas are readily derived from the rate equation from absolute reaction 
rate theory: 


k = Aba exp (- A E/T) = ae exp (4 st /p) exp (-AA H'/RT) 


h h 


The data and the results of the calculations are given in Table VI. 


Table VI: Data and Results Relative to Rate 
of Conformer Interconversion 


° ly ~ lie 
T(°C) (sec) eK) log, , (=F) 
Zo 2.14x10°° eerie: Biome) GO lie) 
20 4.2810” 3.41x10°° Sao GrO4 Oo y=) 
a =a 
0 ro OehG Solos 0 Tao (60,077 =0 506) 
-§8 - 
zleS oo ek) Saisie CLG : 7.72 (+0.14,-0.11) 


A FP" =5.86 + 0.16 kcal/mole 
AH =7.21 + 1.26 kcal/mole 
As" =4.43 + 4.71 cal/mole °K 


4: 
The value calculated for /AF agrees very closely with the value of 


5.7 +0.2 kcal/mole by Windle, Kuhnle and Beck [2]. 
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Figure 9. EPR Spectra — Applied Field Along the nr4 


(eerie i pale xis. “oor Experimcmued 
Bovvom:e simulated. 
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Figure 10. EPR Spectra - Applied Field Along the n+? 


Verrenrct pie ears. Lop: NxXperimemuanw 
Bottom: Simulated. 
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EPR opge tra - Applied Field Midway Between 
the N 2P Orbital and the N-O Bond in the 
eA Prancapal xis Plane. Lone eiroe omental 
BOVGOM: cul aved , 
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Figure 13. EPR Spectra - Applied Field Midway Between the 
X and the Z Principal Axes in the XZ Plane 
(Crystal Axis Syotem Cocrdiva .es+sc=7 775" 
fe=en5).— Mon. 2o-C> Boryon- —lo. ee 
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Figure 14. Arrhenius Rate Plot for Conformer 
Interconversion. 
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Vo Diseusslon 


A. SUCCESS AND Are riICARI Riso Tiberi 

Although the complex model used by program ESRPLT results in 
a very precise prediction of the individual line energies and intensities 
which compose an EPR spectrum, the more approximate approach of 
program NITROX gives plotted spectra which are indistinguishable from 
those produced by ESRPLT. This observation, combined with the fact 
that NITROX can be used to simulate a much wider variety of EPR spectra, 
makes NITROX the better choice as a tool for use in the study of para- 
magnetic system as typified by morpholine N-oxide. In fact, the extent 
to which the simulated spectra are reproductions of the experimentally 
observed spectra could only be enhanced by additional research directed 
toward the determination of separate direction cosine matrices for the 
axial and equatorial proton coupling tensors. 

The possible applications of this experimental tool have been only 
partially explored during the research reported in this paper. Itis 
adaptable to systems containing as many as two or as few as zero 
nitrogen nuclei, and it can include up to six equivalent or inequivalent 
protons in the treatment. It is particularly appropriate, in combination 
with program GEX, to the investigation of exchange phenomena. 

It was observed that simulated spectra are very useful in the 


determination of coupling constants which are not otherwise measurable 
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because of poor resolution. A normal situation was that the i and 
axial proton coupling constants could be measured fairly accurately 
from experimental spectra, but the equatorial proton coupling constants 
usually had to be estimated. A few trials with program NITROX, using 
the measured and estimated coupling constants, usually resulted in the 
specification of a range of uncertainty for the estimated equatorial 


coupling constant of about one gauss. 


Be REGARDING i HYPERFINE COUPLING 

The observation that the average of the principal values of the 
Nik. hyperfine coupling is the same at ~150° C as at 20° C is very 
significant. It implies that the isotropic contribution (Fermi Contact) 
to the A-tensor does not change as the rate of exchange is decreased 
by a factor of at least ten (see section IV: E). In addition, the data 
show that the dipolar interaction increases along the Z principal axis 


and decreases along the X and Y principal axes as temperature is lowered. 


S. SUGGESTIONS FOR FURTHER WORK 

There is still much to be learned about the magnitudes and 
directions of the interactions which give rise to the observed EPR 
spectra of complex paramagnetic systems. The tools developed in this 
research can be used to extract much of this information. Some ideas 
for further research are: 

Ihe Determine the direction cosines for the axial and equatorial 


proton coupling tensors at low temperatures (low rate of exchange). Use 
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these to determine the principal values of the axial and equatorial 
proton coupling tensors. It is expected that this information will be a 
very useful aid to the understanding of the Fermi Contact and dipolar 
interactions. 

De Combine programs NITROX and GEX so that, given the 
energies and intensities of the lines of the base spectrum for the case 
of no exchange, and given the predicted energies and intensities which 
Simulate some experimentally observed spectrum, GEX computes the 
rate of exchange which must occur to cause the base spectrum to evolve 
into the observed one. This should result in an accurate prediction of 
the free energy, enthalpy, and entropy of activation for the process of 
conformer interconversion. 


oe Progress to similar, complex paramagnetic systems. 


D. SUMMARY AND CONCLUSIONS 
i A Useful Set of Tools 
A set of research tools has been developed in the form of 
programs ESRPLT, NITROX and GEX, which should prove very useful in 
the exploration of complex paramagnetic systems. An especially useful 
feature of the spectra produced by NITROX and ESRPLT is that they are 
scaled so as to be superimposable on spectra recorded using the Varian 


SeOCCtIrOMeLer . 
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Zn Comparison of Low and Room Temperature Hyperfine 
Coupling Constants 


Table VII compares observed hyperfine coupling constants 


at room temperature with those observed at -150°C. 


Table VII: Room Temperature versus -150°C 
Hyperfine Coupling Constants for Morpholine N-oxide 


Nucleus Coupling Constant Room Temp. -150°C 
A 8.6G 7.17 G 
14 

N Avy 9.4 G 7 <0 7G 
Aog 348.3 33.0 Ge 

Ane Poo G 23.0 2E€ 

lee? 

axial H AL 11.3G 21.8 G 
Aog ae 23.3 SG 

A 13e6nE (a © 

i 

equatorial H 2 ey 11.3G 22a 
Ao Lie AG 2.3 G 


Notes: 
(1) Hydrogen coupling constants at -150°C are those observed along 

V4 es ; 
the N principal axes. The near isotropy of the room temperature 
hydrogen coupling constants makes the comparison more significant than 
it might otherwise be. 
(2) Exchange at room temperature results in the observation of an 


average value for the axial and equatorial hydrogen coupling constants. 
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Sis Regarding Exchange 


With the aid of program GEX, it was possible to simulate 
Spectrum changes due to exchange phenomena, for spectra recorded in 
the temperature range from -150° to 29°C. The resulting data were 
used to compute the free energy of activation for the process of intercon- 
version, and the data also led to reasonable estimates for the enthalpy 
mucientropy Of activation. The close agreement belween the iree energy 
value obtained in this research and the value obtained from EPR studies 
of morpholine N-oxide in methylene chloride leads to the conclusion 
that the mechanism for exchange is the same in solution as in the 
crystal matrix. 

2 ce Anisotropy 

It was possible, in this research, to observe the effects 
of g-tensor anisotropy in experimental spectra. An example of the 
disappearance of g-anisotropy with increasing rate of exchange is 
detectable in Figure 13. Note the slight asymmetry of the peaks at the 
Opposite ends of the spectrum at -15°C. This asymmetry, which is 
attributable to g-anisotropy, had disappeared by the time the temperature 
nad reached 29°C, as evidenced by the upper spectrum. 

at Alpha 

The magnitude of the value determined for alpha means that 

the angular displacement of the N-0 bond is only five degrees above or 


below some mean plane described by the "average" direction of the N-0 
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bond and the two B carbon atoms. This would imply that only a small 
amount of motion occurs on the N-0O end of the molecule when conformation 


interconversion occurs. 
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APPENDI. 


MECHANICS OF ESRPLT 


Most of the input parameters, and the mechanics of this program 
are covered by comments in the program itself. Several items, however, 
require additional clarification. 

i The Rotation Angles, Theta, Phi and Psi 

The spatial relationship between crystal and principal axes 
may be expressed in a form other than the direction cosine matrix. The 
transformation of one axis system into another can be viewed as con- 
secutive clockwise rotations about the X, Y, and Z axes of the reference 
system by the rotation angles theta, phi, and psi respectively, such 
that after the rotations have been accomplished, the two axis systems 
are superimposed. The rotations about the orthogonal axes of the 
reference system can be effected by a series of matrix operations. 

Once the proper rotation angles have been determined, the transformation 
matrices can be applied to vectors. Consider, for example, a unit 

vector expressed in terms of its components in the crystal axis system. 

If this unit vector is operated on by the proper transformation matrices, 
the resulting matrix will give the coordinates in the principal axis system, 
of each of the components of the vector. Shown below are the proper 
transformation matrices, operating on a vector with coordinates in the 


Poona Zasvsrem Of | 7 d=) 1. 
x Yy Z 


oS 





cosY sin’ 0. cos @ 0 -sin g@ i 0 0 


-sin’ cos¥ 0 De 0 ] 0 aneos. oO sin 8 O 
0 0 l sinf@ OO cosg@ -sin6 cos@ 0 
rotation by ¥ rotation by @ rotation by @ 
about Z axis about Y axis about X axis 

1 1 1 
_* Se _2 
l 0 0 xe x! ; 
x X Y Z, 
xe 0 ] 0 = Gg we a 
y x y Z 
0 0 ] Z' Le Z 
Zz x y Z, 
the vector in the the vector in the 
X,Y,24 system XY Zeiss tem 


ibaat is, ee Ve! a are the coordinates of l in the prime system. 

The three transformation matrices were combined into a 
single, total transformation matrix (TCOS). The rotation angles theta, 
phi, and psi which were used in both ESRPLT and NITROX were determined 
by setting each element of the total transformation matrix equal to the 
corresponding element of the appropriate direction cosine matrix. A 
computer program, LSTPP [9], was designed to calculate 8, @, and ¥ , 
using the method of least squares to give best fit values. 

a The Rotation Matrix RCOS 

A matrix was developed which would simulate the effect of 
an angular displacement of the N-0 bond. This matrix, called RCOS, 
operated on the matrix TCOS. While TCOS was used to express the 


coordinates of a vector in the principal axis system, RCOS transformed 
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the coordinates so that they described the vector in a new principal 
axis system which has been displaced alpha degrees from the old 
one [10]. 
oe PTOP ~ Line Width 
The peak-to-peak line width could be varied as necessary 
depending on the type of spectrum desired. It was found that PTOP = 
0.015 gave good "stick" spectra, while PTOP = 0.25 seemed to work 
best for spectra simulation at -150°C. (Here PTOP is expressed in 
inches. The conversion factor to gauss is 0.06 inches/gauss.) 
ae Scaling 
The spectra which are plotted by this program are scaled 
at 15 inches/250 gauss, which is compatible with the spectrometer 


msead in this research. 
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APPENDIX B 
MECHANICS OF NITROX 

The input parameters to program NITROX are nearly identical to 
those of ESRPLT. The most significant difference is in the meaning of 
imesoarameter SIC (for IC). Im this orogram it gives the Case number for 
the i-th nucleus. The case number specifies the nuclear spin of the 
i-th nucleus, and it also designates the number of equivalent nuclei 
@mine i-th kind. 

ReOS, (COS, PLOP, and scaling areas in EokE Li. The comments 
at the beginning and throughout the program give an adequate description 


of what the program does. 
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